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1 INTRODUCTION 


1.1 BACKGROUND 

This study of voice bandwidth compression techniques was motivated by 
anticipated link margin difficulties in the Shuttle S-band communication 
system. It was felt that by reducing the data rate on each voice channel 
from the baseline 24 (or 32) Kbps to 8 Kbps, additional margin could be 
obtained. Thus, this study was undertaken to determine the feasibility of 
such an alternate voice transmission system. Several factors of prime 
importance that were addressed are: 

1) Achieving high quality voice at 8 Kbps, 

2) performance in the presence of the anticipated shuttle cabin 
environmental noise, 

3) performance in the presence of the anticipated channel error 
statistics, 

4) minimal increase in size, weiaht, and power over the current 
baseline voice processor. 

1.2 TASKS 

The following is a summary of the tasks performed under this contract: 

Task 1: System Requirement Analysis 

Descriptions of the operational environment includinq cabin background 
noise and channel error rates were determined throuah consultation with 
NASA-JSC personnel. Material for system testing was supplied by NASA/JSC. 

It was decided, based on previous work at TRW, to examine various forms of 
adaptive predictive coding (APC) for the Orbiter application. 

Task 2: Compression Technigue Design and Evaluation 

Preliminary analysis of several existing APC algorithms led to the 
selection of a computationally simple algorithm with a noise squelching function. 
This basic algorithm which was developed under a TRW I RAD oroqram was then 
subjected to a preliminary parameter optimization specific to Shuttle require- 
ments. The system was simulated in integer arithmetic on an Interdata 85 
computer, and voice tapes provided by NASA/ JSC were processed through the 
simulation. The properties of the proposed channel were also simulated and several 
tapes were processed throuah the APC system and simulated channel. 
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TASK 3: Parameter Optimization 

Tradeoffs were performed varying the parameters of the chosen system, 
such as frame rate, sampling rate, number of coefficients, etc., and it was 
determined that the original system was near optimum in light of the imple- 
mentation constraints. A revised baseline was established and used to 
process one tape. In consultation with NASA/JSC, It was decided to examine 
a special-purpose hardware implementation based on a micro-controller and 
arithmetic unit. A preliminary design and a sizing were then performed 
for this confiauration. 

1.3 SUMMARY 

An 8 kbps Adaptive Predictive Coding system for potential orbiter use 
was designed and then simulated in integer arithmetic. The system performs 
well and shows good resistance to both channel errors and background noise 
similar to those anticipated in the orbiter application. Channel error 
rates ranging from 10" 4 to 10 -2 were simulated with the result that rates of 
10" 3 or less were judged to have negligible impact on the received voice 
quality, and a rate of 1 0“ 2 , while noticeable, produced no unexpected distor- 
tions. Due to the inclusion of an adaptive squelching unit, the level of 
noise in the received voice wa< often below the input noise level at the 
transmi tter. 

A preliminary estimate of the flight hardware configuration indicates 
that one full duplex system could be implemented in fewer than 200 C-MOS IC's 
(including both LSI and MS I chips) using less than 20 watts. The amount of 
hardware required is insensitive to small variations (ten percent or less) 
in both the data rate and frame size. Thus, the system is easily adaptable 
to minor changes in the Network Signal Processor operation. 
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2 VOICE ALGORITHM 


2.1 THEORY 

The chosen system is a form of adaptive delta modulation (Figure 2.1) 
in which both the quantizer steD size and the Dredictor coefficients are 
optimized over a time interval (called a frame) for the current speech 
statistics. Experimental observations of speech statistics indicate that the 
speech signal can be considered stationary over a time interval of 10 to 30 
milliseconds. Thus, a typical frame time is chosen in this ranae. 

2.1.1 T ransmi tter 

The predictor coefficients are chosen to minimize the power in the 
error, e n> given by 
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n 



M 
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( 2 . 1 ) 


Thus, the power for an N sample frame is 
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Several schools of thought exist on the solution of this system. We 
have chosen the autocorrelation formulation, which yields the simplest solution. 
For this approach it is assumed that S p » 0 for n < 1 and n > N. Thus, 


'S 


HilJ-i I 


5„ V|j- 1 | - "|j. ir each frame 


Is considered independently from its neighbors. We may re-phrase the 
problem as follows: 


(1) Compute 


l i 


N-1 


S n+i 


(2) Solve 


M 

^ aj R|j.,| ■ -R, 1 ■ 


Clearly, the correlation coefficients R ■ form a matrix [r^] which is 
symmetric, has all positive entries, and in which all the elements along the 
diagonal or any off-diagonal are equal. Further, the values R are all 

chosen from the set {R q , . . )• Systems of equations involving such a 
matrix (which is called Toeplitz) are easily solved using a technique known 
as Levinson's recursion^. A flowchart for Levinson's recursion is given in 
Figure 2.2. 

Unfortunately, the solution of this system of equation? may lead to a 
formulation of a filter which is unstable, that is, some of the zeros of 


M 



or equivalently, some of the poles of 
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Finure 2.2 Levinson Recursion 














may lie outside the unit circle. This is most likely to occur due to 
round off and truncation errors when the input speech has one or more 
components with narrow bandwidths. Thus, it is usual to window the frame of 
speech samples with a function which tapers to zero near the end points of 
the frame. A nearly equivalent approach which is computationally simpler, is 
to window the autocorrelation function in such a way as to provide a slight 
increase in bandwidth in the resultino filter. To a qood approximation : If 

a signal S(t) has bandwidth B and autocorrelation function R ( t ) , then the auto- 
correlation function e" alr R( t ) corresponds to a siqnal with bandwidth B + 2 t. 
The appropriate modification to diqital autocorrelation coefficients is 

= exp { -Tt FTi } R i 


where T is the samplinq interval in seconds, and F is the increase in bandwidth 
in Hz. Values of F in the range 10 to 50 Hz have been found experimentally 
to be sufficient to ensure stability. The Levinson recursion is then applied 
to the va 1 ues { R i } . 

A by-product of the Levinson recursion is a set of M coefficients 
M 

(K.) , known as the reflection coefficients. The rronertles of the re- 

J j=l 

flection coefficients have been studied at length. They are bounded by -1 for 
stable filters, and thus provide a stability check. Because of this and the 
fact that a simple algorithm exists for transforming reflection coeff i cieints 
to predictor coefficients, it is clear that it is the reflection coefficients 
which should be quantized for transmission. However, it is important that, the 
predictor coefficients used in both the transmitter and receiver be identical 
in the absence of channel errors. Thus, at the conclusion of the Levinson 
recursion, the reflection coefficients are quantized, and the quantized 
values are applied to a short form Levinson procedure (Figure 2.3) to produce 
a set of predictor coefficients for use in the transmitter's quantizer loop. 

Another by-product of the Levinson recursion is the quantity, a, which 
measures the power of the prediction error. The rms prediction error per 
sample is <e n > = /a/N. This prediction error qives a good indication of 

the minimum quantization error |)2>l/2. 
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Figure 2.3 Conversion of Reflection Coefficients to Predictor Coefficients 







With reference to Figure 2.1, it can be shown thct the optimum 
quantizer setting, Q, which minimizes quantization error, is given by 

o -<ig> 

For most voice signals, it has been observed that for reasonable choices 
of quantizer level , Q, 

< i *„ i > *- 7 <^>' /2 

As the system tracks the J nput siqnal, S n , the siqnal r n is forced to follow 
the input signal. Thus, che prediction error, 6 , within the loop approxi- 
mates the theoretical prediction error, e n . So a qood approximation to the 
quantizer level is 


Q ■ .7 <*> 1/2 * - 7 ^ 

It has been found experimentally that the ratio a/R of the output power 
of the ideal predictor to the input power does not exceed 0.36 for voice 
signals. Thus, this ratio is also used in computing the quantizer level, i.e., 

Q = .7 F( a /R 0 ) 

when a/R o approaches 1 indicating severe noise in the input, F(a/R Q ) approaches 
zero, thus reducinq the ability of the loop to track the siqnal. 

In the loop, the output of the quantizer is fed back through the predictor 
thus providing a closed loop prediction system. In addition, the quantizer 
outpjt (values of -Q only) is converted to a bit stream called the residue 
at a rate of one bit per sample interval. The residue for each frame is 
multiplexed with the corresponding coded parameters to aive the channel bit 
stream. 
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2.1.2 Receiver 


In the receiver or demodulator, the residue bit stream is converted 
to a stream of -Q values and used to drive a synthesis loop (see Figure 2.1). 
The predictor coefficients which are fixed "ver a frame, are obtained from the 
received reflection coefficients by the short Levinson recursion (Figure 2.3). 

2.1.3 Parameter Encoding 

The quantizer level has a wide dynamic ranqe. Thus, it is coded log- 
arithmically. Subjective listeninq indicated neol'qible difference between 
quantization to 4 bits and no quantization. This 4-bit coding scheme was 
used for all simulations and testing. 

The reflection coefficients were coded linearly to 4 bits each. Aqain, 
subjective listening tests indicated little if any difference between this 
quantization scheme and no quantization. Other quantization techniques, such 
as log bilinear, have been shown to yield a slightly hiqher quality voice in 
other types of compression systems. It was felt, however, tnat the siiqht 
improvement that miqht be obtained here did not justify the additional 
computational complexity. 

2.1.4 Baseline System 

Preliminary subjective listeninq tests we^e used to select the followinq 
baseline system: 


Sampling Rate: 

7 KHz (12 1 

bi ts/sample) 

Frame Time: 

20 ms 


Parameters: 




Q Coded 

loqari thmical ly to 4 bits 


K 1 Coded 

linearly to 4 bits 


<2 Coded 

linearly to 4 bits 


Coded 

linearly to 4 bits 


K^ Coded 

linearly to 4 bits 



Tota 1 : 20 bi ts/frame 

Data Rate: 

Parameters 

1 Kbps 


Residue 

7 Kbps 


Total 

8 Kbps 
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This is the system that was used for simulation and testing. Detailed 
tradeoffs of these parameter values are discussed in Section 2.4. 

In practice, several bits per frame could be devoted to frame ryr h. 

A typical allocation with provision for frame synch is: 

Sampling Rate: 6750 

Frame Time: 20 ms 

Parameters: 20 bits/frame 

Synch: 5 bits/frame 

2.1.5 Channel Characteristics 

In addition to resistance to background noise in the speaker's environment, 
the resistance of the system to channel jrrors is of importance. In ord^r to 
determine the effects of the proposed channel, the statistics of burst errors 
re: •• . <ng those of a rate 1/3 Viterbi decoder with constraint length 7 
yi*. * .j 'jO veloped for error rates of 10 10 and 10 . Appendix A describes 

tne channel simulation techniques in detail. 

2.2 INTEGER SIMULATION OF NASA-APC ALGORITHMS ON INTERDATA 85 

An integer simulation of the NASA-APC algorithms was developed on the 
Interdata 85, a 16-bit machine, in two versions: 

IAPCDC: Disk-to-disk 1/0 

IAPCTC: Tape-to-tape 1/0. 

The programs are identical except for the 1/0 subroutines called. 

The programs execute NASA-APC transmitter and receiver algorithms, and 
simulate burst errors (as specified by parameter input) on the communi- 
cation channel . 

The main programs are in FORTRAN and "se functions and subroutines to 
simulate fractional integer arithmetic. The simulation equated + 1 to 32767 
and expressed fractions as a proportion of the + 1 base. 

The programs require the following parameters during initialization: 
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M 

Number of Coefficients 


SR 

Sampling Rate 


NSPF 

Number of Samples per Frame 


ERRATE 

Channel Error Flaq 

+, Simulate channel errors 

-, Do not simulate channel errors 


FRAC 

Quantization Level Gain Control 


F 

Quantization Level Noise Squelch Control 


JAQ. 

Reflection Coefficient Quantization, i * 

1, M 

JREC 

Number of Input Frames to Skip 3efore Process 

NUMREC 

Number of Frames to Process 
-, Process to enH-of-file 


NPC 

Maximuti. Channel Burst Length 


EC 

Channel : E s /N 0 


PC i 

Probability of Burst of Length 1,1*1, 

NPC 

IT i 

Input Seeds to Ranuom Number Generator 
i = 1 , 19 



At the end of initialization, the programs inform the operator of the 
bit rate. If channel errors are simulated, the program prints the total number 
of errors simulated every 100 frames. When the number of frames designated 
h;.s been processed, the proqram prints the number of uses of each burst lenqth 
and the calculated error rate. A flowchart of the program is given 
in Figure 2.4. 
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FIGURE 2.4 INTEGER SIMULATION OF APC 








2.3 EXPERIMENTAL RESULTS 

Although various signal-to-noise ratio measurements were made on the 
chosen system (see Section 2.4), .1 much more meaningful measurement of 
performance was obtained through standardized intelligibility tests per- 
formed at Fort Huachuca. These tests were used to determine the intel- 
ligibility in the presence of channel errors as well as the intelligibility 
over an error free channel. The channel simulation used is discussed in 

Appendix A. Figure 2.5 shows the results obtained for the four channels. 

.3 

As can be seen, error rates as high as 10 have virtually no effect on 
intell igibil ity . 
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2.4 TRADEOFFS 


In order to optimize the chosen system, several parameters were 
varied and the signal-to-noise ratio in the transmitter loop was measured 
on a frame basis. The SNR measurement is defined as 

SNR * 10 loq ln [ Input Signal Power 

lu |_ Quanti za tion Noise Power 

The quantization noise power and the input siqnal power were each avsraqed 
over a frame, and then the SNR for that frame was computed as above. The 
parameters examined were frame length, number of coefficients, bit alloca- 
tion (quantization) of coefficients, and the quantizer level scale factor. 

The performance was evaluted by examininq the signal -to-quanti zer 
noise ratios on 0.6 seconds of speech. The input SNR of the diqitized 
speech was estimated at 30 to 40 dB. 

Three frame sizes were examined and the -esults are shown in Figure 2.6 
For the most part, the shortest frame (15 ms) qave sliqhtly hiqher SNRs on most 
of the samples. In one reqion of very rapid transition, the 20 ms frame 
gave better results than either the 15 or 25 ms frames. This was probably 
due to the particular alignment of each of the frame boundaries at that instant. 
Listening indicated no perceptual difference between the three variants. 


Figure 2.7 shows the results of comparinq three systems: one with 3 

coefficients, one with 4, and one with 6. Both the SNRs and listeninq in- 
dicated sliqht degradation with the use of 3 coefficients, and virtually no 
difference between the systems with 4 and 6 coefficients. In all three cases 
each coefficient was quantized to 5 bits. 

I 

Since little is gained by using more than 4 coefficients, we next examined 
bit allocations for 4 coefficients. A total of 16 bits were allocated in four 
different ways: 


Allocation 


1 

2 

3 

4 


4 

5 
5 
5 


4 

4 

5 
5 


4 

4 

4 

3 


4 

3 

2 

3 
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Figure 2.7 


The first allocation was judqed to be slightly worse than the remaininq 
three which were all roughly equivalent. Figure 2.8 gives the results 
for the first and fourth allocations. 

The final parameter examined was a scale factor used in setting the 
quantizer level Q, i.e., the value of Q before quantization was 


F (a/R 0 ) 


F(a/R) was taken as (1 - .9 5 —). Thus, the theoretically optimum value of f 

K o 

is 0.73 for errors with Laplacian statistics and 0.82 for errors with Gaussian 

statistics. In Figure 2.9 the resulting SNRs are shown for four values of 

f ranginq from 0.5 to 1.1. Although the SNRs indicate that the values 0.7 

through 1.1 are roughly equivalent, and 0.5 is superior only on near silence 

intervals, in fact 0.5 was judqed perceptually the most pleasing on all speech. 

T31 

This is probably rit.e to a phenomenon which has been noted elsewhere 1 - , i.e., 
that the human ear tolerates the distortion due to slope overload more readily 
than quantizing noise. In further perceptual evaluation, values of f in the 
range 0.4 to 0.5 were found to be the most qenerally satisfactory. 

As a result of the above evaluations, the following alternate baseline was 
establ i shed: 


Sampling Rate: 7.0 KSPS 

Frame Length: 25 ms 

Parameter Encoding: 


Total 

Synch: 


l <2 5 bits 

K 3 3 bits 

K 4 3 bits 


20 bi ts 
5 bits 
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Parameter Rate: 
Synch Rate: 
Residue Rate: 


800 bps 
200 bps 
7000 bps 

Total Data Rate: 8000 bps 

One tape, an FM news broadcast, was processed using both this system 
and the baseline used for the other tests. No significant perceptual 
differences were noticed between the two systems. 

The conclusion we reach is that there is substantial flexibility i r. the 
allocation of bandwidth in the neighborhood of the original baseline alloca- 
tion. Minor changes resulting, for example, from the Inclusion of either 
synch or error protection a few parameters would not significantly influence 
received voice quality. 
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3. IMPLEMENTATION 


3.1 NETWORK SIGNAL PROCESSOR INTERFACE 

The two voice compression modems utilized in the orbiter NSP are shown 
in Figure 3.1. The analysis and synthesis functions of eacii modem are 
included in a single package, called a 1-AU, to allow sharing of common 
processing elements. Within the NSP, the voice compression units interface 
directly with the TDM MUX/DEMUX equipment where the compressed voice is 
exchanged. Timing information for external framing and clocking the digitized 
voice into and out of the voice units is provided by the MUX/DEMUX units. 

Within the units frame synch for the voice algorithm can be provided, or the 
NSP may provide frame synch signals. The choice between these two options 
should be based on a study of the overall impact on the NSP hardware (see also 
Section 3.3). Each unit interfaces with the audio center via two analog voice 
signals: one from ground to orbiter, and one from orbiter to ground. The data 

shown are for 8 kbps APC, or the 32/24 kbps VSD currently baselined, which is 
shown parenthetically. 


\ 



Figure 3.1 Voice Modems and Network Processor 
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3.2 NASA-APC IMPLEMENTATION ON THE 1-AU PROCESSOR 


A full-duplex NASA-APC system has been designed for real-time imple- 
mentation utilizing the 1-AU processor shown in Figure 3.2. A simplified 
block diaqram of a hook-up between two voice processing systems is Illustrated 
in Figure 3.3. The system is desiqned to start randomly when powered up 
and run continuously until powered down. 

A preliminary coding exercise was performed to provide detailed timing 
estimates and hardware requirements. The timing estimates were developed 
on a 300 ns instruction execution time base, and are summarized in Table 3.1. 

The 1-AU processor does not have normalize or divide instructions. The 
APC transmitter algorithms require 1 normal ization and 5 divisions per frame. 

It was determined that these operations could be implemented in the existing 
software capability without additional hardware. 

The APC algorithms used in the integer simulation on the Interdata 85 
were converted to the 1-AU processor in a real-time frume-work augmented by 
the input/output operations necessary to the system. The sampling rate was 
modified from 7030 samples per second to 6750 samples per second to accommodate 
both APC data and synch in the 8000 bits per second communication channel, since 
this is a worst case configuration. The flowchart presented in Figure 3.4 
provides a functional summary of t he program. 

The input/output of the NASA-APC transmitter is asynchronous to that of 
the receiver; however, the proqram will function synchronously. Current 
transmitter output and receiver input are buffered in the MUX buffers, while 
the next frames are being processed. The D/A outputs are double-buffered 
in the RAMs as are the A/D inputs. The latter are double-buffered in both 
RAMs to provide maximum efficiency durinq the construction of the autocorrelation 
matrix (i.e., both the multiplier and the multiplicand can be loaded with 
one instruction). 

Notification of input/output status is via discretes. The execution 
of the APC algorithms is time-sliced accordingly to ensure proper servicing 
of the A/D and D/A sampling rate. At a sampling rate of 6750 samples per second, 
the program is structured to check the 1/0 status of the A/D end D/A within 
148 micro time increments. In the event the A/D status indicates no input 
available, zeros will be processed. 
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Table 3.1 


NASA-APC Implementation on 1-AU Processor 
Timing and Sizing Estimates 


Function 

Op. Cycles 

RAM 

ROM 

Construct Autocorrelation Matrix 
and Apply Stability Weighting. 

3553 

146 

123 

Invert Autocorrelation Matrix 

372 

7 

132 

Calculate Ouantizer Level 

136 

22 

66 

Quantize Reflection Coefficients 

16 

4 

( 

Convert Reflection Coefficients 
to Prediction. 

26 

4 

26 

Execute Predictor Filter 

3809 

145 

229 

Pack Transmitter Output Buffer 

945 

10 

34 

Unpack Receiver Input Buffer 

945 

10 

34 

Determine Quantizer & Predictor 
Coefficients . 

30 

- 

30 

Reconstruct Signals 

2185 

145 

85 

Clip D/A Output 

407 

- 

6 

Refresh RAM Constants 

360 

- 

6 

A/D, D/A Input/Output 

540 

540 

540 

Service MUX Buffers 

40 

\ 

10 

Housekejp RAM Buffers 

280 

- 

10 

TOTAL 

13644 

1033 

1347 
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Figure 3. 





[CALCULATE QUANTIZER LEVEL j 

[ QUANTIZE REFLECTION COEFFICIENTS! 

CONVERT REFLECTION COEFFICIENTS 
j TO PREDIC TION COEFFI CIENTS | 


EXECUTE PREDICTOR FILTER 


PACK TRANSMITTER OUTPUT BUFFER - j 


UNPACK RECEIVER INPUT BUFFER j 
| LOOK-UP QUANTIZER VALUE j 


RECONSTRUCT SIGNAL 


CLIP D/A OUTPUTS | 


REFRESH RAM CONSTANTS 

c 

COMPLETE A/D-D/A INPUT/OUTPUT 


SERVICE MUX BUFFERS 


HGUSEKEEP RAM BUFFERS 


4 Functional Flow Chart of Frame Processing 
of NASA-APC Implemented on 1 AU Processor 
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When the APC algorithms nave been executed, the remaininq frame time is 
spent completing the A/D and D/A input/output, switching buffer references, 
and housekeeping the MUX buffers. If the status of the receive^ *UX buffer 
indicates no data available, zeros will be processed. 
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3.3 FLIGHT HARPWARL CONFIGURATION 


The functional configuration assumed for flight operation is shown in 
Figure 3.5. The framing shown corresponds to the baseline described in 
Sections 2.1.4 and 3.2. The APC analysis and synthesis would be accomplished 
on a 1-AU unit as descr.ued in Section 3.2. 

A preliminary IC count including MSI and LSI is: 


Frame Synch & NSP Interface 20 ICs 
Serial/Parallel 4 Parallel/Serial 6 ICs 
1-AU 158 ICs 
A/D, D/A & Buffers 10 ICs 

Total 194 ICs 


Since the sizing in Section 3.2 indicates 20* bus* •»r,d 80% dead time 
at a 300 ns cycle time, a cycle time of 1 psec would yield a 68% busy status. 
This would permit implementation in C-MOS. Fiqurinq from a 100 mw per chip 
average power dissipation for C-MOS, this gives a total power consumption 
of roughly 20 watts per full duplex voice unit. At standard flight hardware 
packing densities of 55 IC per 20.32 cm x 15.24 cm (8"x6") board, tour 
boards would be required. These four boards may be mounted in a module, 
roughly 5.08 cm x 15.24 cm x 20.32 cm (2"x6"x8") weighing 1.1 kg (2-1/2 lbs.). 

The only effect of removing frame synch, i.e., assuming it is provided 
by the NSP, is to decrease the chip count to 175 ICs. 
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4. RECOMMENDATIONS 


4.1 BASELINE 


The following allocation represents a base point about which variations 
are possible to accommodate synch, error protection, or minor changes in the 
NSP baseline: 

Sampling rate 7 kbps 
Frame rate 50 Hz 
Parameters: 


K t : 5 bits 

K 2 : 5 bits 

K 3 : 3 bits 

IC,: 3 bits 

Q : 4 bits 

RESIDUE: 7 kbps 

PARAMETERS: 1 kbps 

TOTAL DATA RATE: 8 kbps 


4.2 FUTURE WORK 

Tests conducted on similar systems have indicated that the received 
voice quality is particularly sensitive to errors in certain of the parameter 
bits. Thus, an obvious possibility which should be examined is the alloca- 
tion of a portion of the allowed data rate to error correction. This could 
best be accomplished by first determining the relative sensitivity of the 
received voice to errors in each of the parameter bits and then applying 
coding beginning with the most sensitive bits and proceeding toward the least 
sensitive until an optimum is found for a typically noisy channel. Since 
the optimum split between information and redundancy rates is probably 
dependent on channel characteristics, different techniques might be used on 
the uplink and downlink. 

Another area worthy of investigation is the portion of the algorithm 
which sets the quantizer level. Results obtained so far indicate that this 
is the most sensitive part of the algorithm, and that an optimum technique 
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has yet to be developed. A systematic examination of the effects of the 
quantizer level on various signal types is needed to derive a truly optimum 
setting. Once this optimum setting is derived the technique woulr be 
refined to provide noise squelch specific to the orbiter environment. 
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APPENDIX A. SIMULATION OF CHANNEL ERRCR STATISTICS 


It was assumed that the APC bit stream would be transmitted using an 
optimum rate 1/3, constraint length 7 convolutional code. Dr. Gaylord Huth 
of Axiomatix provided an algorithm for generation of upper bounds on the 
probability of burst o' errors as a function of E s /N Q and burst length. The 
probability of a burst of length b can be upper bounded by 

^b, max 


V I 


n.. P. 

bk k 


k=k 

K K b, min 


whe~e P^ is the probability of incorrectly choosing a weight k code word over 
the all zeroes code word, and n^ k is the number of weight k code words 
caused by bursts of length b. P k is given by 


P 


k 





dt 



n^ depends on the code used and there is no general formula for computing it. 
For the code assumed here 


1 j = 14 

n l, j 

0 else 

1 j = 16 

n 2, j = 

0 el se 

The Fortran version of an algorithm supplied by Axiomatix is shown in Figure 
A.l. This program was used to generate values of n^ k for b > 2, and for 
all values of b. 

The total probability of error at the output of the Viterbi decoder can 
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be upper bounded by 


p e " Q i + vH 1 * 2) 

b=2 

This upper bound assumes that a burst begins and ends with a "1" and contains 
a random binary sequence between the two end points. In fact, a burst as 
defined for this code contains no more than six consecutive zeroes at any 
point internal to the burst. Thus, P e as given above is upper bounded in 
both and its weighting factor. An approximation to P g is also computed 
by the program in Figure A.l as 

N 

p e = Qi + y + 2) 

b=2 

where N is the maximum value of b used. 

Due to the excessive execution time of the program for burst lengths 
onger than 20, it was not practical to compute values of beyond that 
point. However, for E $ /N o - 0.0 [i.e., P g a 10’ 3 ], bursts of these longer 
'lengths are nearly as probable as shorter bursts. Efforts to approximate 
the probabilities of longer bursts were not successful, hence the curves for 
versus b were extrapolated to larger values of b until the value of P £ was 
approximately the predicted error rate. Table A.l gives the values of Q^, 

P g , and the probability of error observed in the simulation for E $ /N o = 0.52, 
O.b, and 0.7. Figure A. 2 snows calculated and extrapolated values of 0^. 

The routine for simulating the burst errors acts as a two state device 

where the first state represents the generation of a string of zeroes and the 

second state represents trie generation^of a burst. The burst generation 

state is entered with probability P = ^ Q. where N is the maximum burst 

b-1 D 

length being used. The distance between bursts is a random variable uniformly 
distributed on (0, 2/P). Thu:, the average distance between bursts is 1/P. 
When the burst generation state is entered, a burst length, .1, is randomly 
chosen according to the conditional distribution {Q^/P}. Then a "1 is 
generated followed by J-2 "0"s and l: l's which are equally likely unless 
six "O' s have been generateu in succession, at which time a "1" is forced. 
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Table A.1 


Calculated and Observed Probabilities 



E s' N o 

« .52 

E s /h 

'o 

* .6 

V N o • - 7 
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"b 
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lb 
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«b 
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00 

10“ 4 

.25 

X 
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2 

.23 x 
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.60 
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3 

.30 x 
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.76 
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10-5 

i 

O 
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.18 
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The Jth bit In the burst Is forced to be a "V. After the Jth bit has been 
generated the routine returns to the zero generation state 
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